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Abstract 

Prom a picture of the X(3872) where the resonance is a bound state of DD* — c.c, we evaluate 
the decay width into the J/ip^ channel, which is sensitive to the internal structure of this state. 

For this purpose we evaluate the loops through which the X(3872) decays into its components, 
and the J/ijj and the photon are radiated from these components. We use the local hidden gauge 
approach extrapolated to SU{4) with a particular SU{4) breaking. The radiative decay involves 
anomalous couplings and we obtain acceptable values which are compared to experiments and 
results of other calculations. 

Simultancusly, we evaluate the decay rate for the X(3872) into J/ipoj and J/ipp, and the results 
obtained for the ratio of these decay widths are compatible with the experiment. 

We also show the grossly unacceptable results that come from taking only the D^D*^ + c.c. 
component. 



I. INTRODUCTION 



The first observation of tlie X(3872) decay into J/ip'^ was reported by tlie BELLE col- 
laboration in pjj. Later on this decay mode was confirmed by the BABAR collaboration in 
[2] and more recently again in BELLE in |3|. Theoretically this decay mode had already 
some early attention and was studied in |5-[7j assuming either a charmonium state or a 
molecular state. A thorough discussion of the different models used and the results can 
be seen in [8], which has also been updated recently in [9]. A recent work assuming the 
X(3872) to be a charmonium state is presented in [10] and assuming it to be a tretaquark 
in [11]. In [12] it is assumed to be a mixture of a charmonium and a molecular component, 
and uscing QCD sum rules a good rate is obtained for the J/ip'~i decay mode versus the 
J / ijjTi'^ T{~ one, which is evaluated in [13]. In [8] the authors consider, like in the work of 
[5J, the X(3872) resonance to be a molecule of D^D*^ — c.c. and, in addition, they include 
the possibility of a cc admixture. In [S] an effective Lagrangian is postulated to provide the 
coupling of the X(3872) to the D^D*^ components, with an unknown wave function. The 
effective coupling needed in the loops for radiative decay of the X(3872) is obtained using 
the Weinberg compositness condition [HI [IS], reformulated in [16] as = — (^G)"^, where 
G is the loop function of the and Z)*° propagators. The procedure has been shown to 
provide a fair description of the molecular states in other works [T714T9] . The results of [8] 
are tied to unknowns on the regularization of the loop functions, the Am parameter used in 
[S], and the binding. The results obtained for the X(3872) decay into the J/ip'j channel are 
of about 125-250 KeV taking reasonable values for the Am parameter between 2 and 3 GeV. 

In [9] the authors include the charge components of D^D*~ — c.c, which were found 
necessary to explain the ratio of X(3872) to J /ipp and J/ipu in [161 EO]. The novelty with 
respect to the previous work of [8] is that the authors use a smaller Am cut off, of the order 
of 0.5 GeV to regularize the loop function, such that the wave function of the D^ID*^ — c.c. is 
much more extended in space. The final results of the new evaluations differ quantitatively 
from those of [8] and are now in the range of 2-17 KeV. It is then clear that a more systematic 
approach to the problem has to be done if one wishes to obtain accurate numbers from the 
molecular picture of the X(3872). This is the purpose of the present paper. 

A consistent dynamical picture of the X(3872) in the coupled channels of DD* — c.c. was 
elaborated in [20j| using an extrapolations to SU(4) of chiral Lagrangians used in the study 
of pseudoscalar meson interaction with vector mesons [21]. This is equivalent to extending 
to SU (4) the local hidden gauge approach of [221425] with a particular SU{A) breaking. 
Given the subtlety of the small binding for the neutral D^D*^ — c.c. component versus 
the about 7 MeV binding for the charged D~^D*~ — c.c. components, a coupled channel 
approach considering these explicit channels with their exact mass and not assuming isospin 
symmetry, was done in [26] , concluding that the coupling of the resonance to the neutral and 
charged components was very similar, which tells us that in strong processes the X(3872) 
behaves as a rather good 1=0 object. The D^D*~ — c.c. components were also included 
in [201 [25] . On the other hand, the study of the wave functions of the resonance for each 
of its channels indicated that the couplings provided essentially the wave function at the 
origin, which indicates that the neutral component, extending further away in space than 
the charged components, because of the small binding, and having the same strenght at the 
origin, will have a larger probability in the wave function. But this has not repercussion in 
strong processes, of short range, or electromagnetic ones where there is production of real 
photons radiated from the components of the wave function, because they see the strength 
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of the wave functions at the origin. In sum, the charged D~^D*~ — c.c. components have to 
be considered in the study of these processes. 

In the present work we follow the approach of [151 [2S] where all the couplings are accu- 
rately determined from the unitary coupled channel approach and are tied to the binding of 
the X(3872), which is generated dynamically as a composite state of DD* in this picture. 
The mechanisms for radiative decay are then basically the same as in [9], except that we 
also have contribution from the DgD* components and have, although not much, different 
couplings of the resonance to the neutral and charged DD* components. The work is also 
technically different. Our approach has not ambiguities about the regularization of the loops, 
and most of the terms are shown to be convergent. Some terms are formally divergent, but 
we can isolate the divergence into a term proportional to the same loop function G which 
appears in scattering. The function G is regularized in the scattering problem in order to fit 
the position of the resonance, so when it comes to evaluate the radiative decay it is already 
fixed. This makes the scheme fully selfconsistent and predictive, since one does not have 
to rely upon unknown parameters that have proved to have a strong repercussion in the 
numerical results from the former studies. 

Traditionally the X(3872) could be considered as a = 1^+ or = 2^+ state, and 
there is a work similar to the one of [9] but assuming = 2 ^ [27] . Here we will continue 
to use the = which is supported by recent analysis of data in [281 129] . 

Our work proceeds as follows: in the next section we present the formalims for the 



work with the Feynman diagrams used and the scheme to evaluate them. In section III 



we present the results for X(3872) — )■ J/ip'y, J/ipoo, J/ipp and compare them to experiment. 



discussing the role of the charged components of the X(3872) wave function. In section IV 
we summarize our results. 



II. FORMALISM 

In [201 EQ], the interaction between pseudoscalars and vector mesons is studied includ- 
ing the charm sector. The potential is like the Weinberg-Tomozawa interaction between 
pseudoscalar mesons but including the vector meson fields [21] . In [201 126] , all the different 
currents within the SU{4:) scheme are classified in terms of SU{3) currents, and the break- 
ing symmetry parameters are introduced to account for the suppression of the heavy meson 
exchange. Within this formalism, the X(3872) is a dynamically generated resonance from 
the interaction of DD*, having an eigenstate of positive C-parity with isospin / = 0. It also 
has some component of DgD*. In fact, the basis of positive C-parity and / = for these 
two channels corresponds to: 

-^\{D*D - D*D),I = 0, J3 = 0) = -\{D*+D- - D*- D+ + D*^ D^ - D*^ D^)) 

f \ m 

- D:D,),I = 0,/3 = 0) = -l^VD'^D- - D'-Dt)} . 

In [2n] it was found that the X(3872) had couplings to the charged and neutral components 
of DD* that were very close to each other, implying an approximate / = character for the 
state. Since the masses and bindings used in [201 ES] have been updated, we have redone 
the calculation of [201 ES] with updated masses, assuming the present binding of 0.2 MeV 
of the X(3872) with respect to the D^D*^ — c.c component. The result of the couplings are 
shown in table HI 
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Channel 


\9R- 


,pv\ [MeV] 


{K-K*+ 


-c.c.)/V2 


-53 


{K^K*^ - 


-c.c.)/^ 


-49 


{D-D*+ 


-c.c.)/V2 


3638 


{D^D*^ - 


-c.c.)/V2 


3663 


{D-DI+ 


-c.c.)/\/2 


3395 



TABLE I: Couplings qr of the pole at (3871.6 - iO.OOl) MeV to the channels {an = -1.27 

here). Table taken from Ref. [2^] . 

From the couplings in Table [T| we observe that there is some isospin violation, which is 
however very small, less than 1%. Intuitively, one might think that the D^D*^ component is 
the only relevant, because the binding of the DD*^ is very small, of the order of 0.2 MeV and 
the wave function extends much further than for the charged component, which is bound by 
about 8 MeV. However, as we mentioned, the relevant interactions in most processes are 
short ranged and then the wave functions at the origin, proportional to the couplings, are 
what matters. In this sense it is found in [201 [26] that in the limit of the D^D*^ binding going 
to zero, all couplings go to zero, but the ratio of the couplings of the charged and neutral 
components goes to a constant close to one, which guarantees that the charged component 
will play an important role in physical processes. The conclusion is that the wave function of 
the X(3872) is very close to the isospin / = combination of D^D*^ — c.c. and D^D*^ — c.c. 
and has a sizable fraction of the D^D*^ — c.c. ofEq. 0. 

From this table we can also see that the couplings to the K~K*^ — c.c. and K^K*^ — 
c.c. channels represent less than the 1% of the contributions from the other channels (the 
TT^p'*" — c.c. has even smaller strenght). Therefore, we will treat the X(3872) as if it were 
dynamically generated only from the last three channels in Table |Tj 

In this picture, the X(3872) decays into J/^'y through the diagrams shown in Fig. [!} 
From this figure we observe that there are four kinds of different Feynman diagrams, all of 
them with an anomalous vertex coupling two vectors and a pseudoscalar (VVP), depending 
on whether the diagram contains a PPV or a 3V vertex, or the photon emerges from the 
anomalous vertex. To begin with, there are three different channels: D^D*^, D^D*^ and 
D^D*~ , which lead to 12, plus another 12 for the complex congujate, Feynman diagrams to 
evaluate. The formalism used is very similar to the one of [3T| where the authors study the 
radiative decay of the dynamically generated resonance K2{14:30) [52] into K'y, via diagrams 
containing anomalous vector-vector-pseudoscalar vertices. The VPP, 3V and V7 vertices are 
evaluated using the local hidden gauge approach [221 - I25] which automatically incorporates 
vector meson dominance, by means of which the photons couple to other hadrons converting 
itself into u, (p and J/ip. As a consequence of this, we are also able to evaluate the rates 
of the X(3872) decay into J/ipp, J/ipco and the ratios of the decay rates, which can be 
compared to existing data. 



4 



Jl^{p) 



X{P) 




\ Pi{q-p) 



ViP-q) 



(1) 



► — ^AA/ 7(fc) 




J/*(p) 



P(g) ^ »^ ^AAy 7(fc) 




» ^AAy lik) 




V{P-q) 



J/*(p) 



(3) (4) 

FIG. 1: Different types of Feynmann diagrams for tlie decay of the X(3872) into J/ipl- 



In summary, the Lagrangians we need in order to evaluate the amplitude are listed below: 

G' 



Cv, = -Mf-A,ly^Q) 
9 

Cppv = -tg{V^[P,d,P]) 



(2) 

(3) 

(4) 
(5) 



with e the electron mass (eV47r = a), G' = Sg'^iAn^f), g' = -GyMp/ {y/2f), Gy = 
f/V2 and g = My/2/. The constant / is the pion decay constant = 93 MeV, Q = 
diag{2, —1, — 1,2)/3 and My is the mass of the vector meson, for which we take Mp. 

The P and V matrices contain the 15-plet of the pseudoscalars and the 15-plet of vectors 
respectively in the physical basis considering rj, rj' mixing |33] : 



P 





7r+ 


K+ 










D- 


K~ 










D+ 


Dt 


Vc 



(6) 
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and represents the vector nonet: 



*0 \ 



V2 



p- ^ K*^ D*- 



K*- K 



*o 



v 



DT 



J/i'J 



(7) 



From Eqs. Q, Q, Q and Q we can write the vertices involved in the diagram of type 
;i) of Fig. [l|as 



tRVP — 9x ^ 



t 



Vp-t 



P M 



e) 'e 



t 



pp,J/4> 



Pv g {2q - p)^t 



(8) 



where gx = 3638/^2, 3663/^2, 3395/^2 MeV, for D-D*+, D^D*^, D*+ and -3638/^2, 
-3663/^2, -3395/^2 MeV, for D+D*-, D^D*^, D+D*^- respectively and P, Py and A are 
numerical factors. 

The — )■ 7 conversion essentially replaces, up to a constant, e^^*''' by e^'^'*. Therefore, we 



can write the amplitude of the diagram (1) depicted in Fig. [T]as 

- iti= -BegxG' j 



1 



1 



X 



q^ — im?p {q — pY ~ "^p {P " ^lY " ' 



(9) 



where B = PAPy (the values of B for each case are shown in Table [ITj). Summing over the 
polarizations of the internal vector, we have 



E 



JV) (V) 



-gpp' 



{P-q)p{P-q) 



(10) 



When contracting with the antisymmetric tensor e^^'^^ , the term [P—q)p[P—q)pi disappears. 
Thus, we have an integral like 



{2q-p)^{p + k- q)a 



(27r)^ (g2 — w?p + ie)((g — pY — rn?p^ + ie){{j) + k — qY — my + ie) 
= i{cigp,a ~\~ bkfj^koi + cpak^ + dkaP^j, + epaP^) (11) 

because of Lorentz covariance. After contracting with the antisymmetric tensor e"^"'^ and 
applying the Lorentz condition p^e^-^^^^'^ = 0, only the coefficients a and c remain to be 
evaluated. The a coefficient is related to the logarithmically divergent part of the integral 
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in Eq. (11) and therefore the evaluation of this coefficient needs a special treatment as we 
will see later on. We arrive to an amplitude of the form 

ti = BeG'gxe''^^\at^^'^^ + cp^k ■ e^"^^)efh,e^^^ . (12) 

Now we want to evaluate the a and c coefficients. We do it using the formula of the 
Feynman parametrization for n = 3, 



a/37 



2 / dx I dy 



1 



^ [a + (/3 - a)a; + (7 - P)y]^ 



(13) 



In the integral of Eq. (11), we can perform the above parametrization with 

a = {q-pf -ml^ 
4 



(3 = — mj 



J = {p + k — q)'^ — m 



V 



(14) 



We define a new variable q' = q + p(x — y — 1) — ky, such that the integral of Eq. ( 11 ) can 
be expressed as 



4: dx dy 



d'^q' {q' +p{l- x + y) + ky)^{k - q' - p{y - x) - ky). 



(2vr) 



iq" + si) 



(15) 



with 



Si = —mp^ + {mp^~mp)x + {k +mp — mY)y+p {x — y){l — x + y) + 2pky{x — y) — k y . (16) 



From Eq. (15), we must take the iag^a and icpak^ terms. The c coefficient can be evaluated 
very easily, since 



d'^q' 



(g'2 + s, 



2ii 



(17) 



and we have 



1 

8^ 



dx I dy^^ — — 
Jo si 



(18) 



The evaluation of the a coefficient is a little bit more elaborated. We have the identity 

dS' 



p1 px 

iagtia = -4 / dx dy 
Jo Jo 



(27r)4 (g'2 j^si + ie) 



(19) 



and after taking the trace 



la 



dx / dy 



J2 







(27r)4 (g'2 + + lef 



(20) 



This part is logarithmically divergent and we will relate it to the two-meson function loop 
G{P): 

d\ 1 1 



G{P = p + k) = i j 



{2t:Y g2 _ _|_ _|_ _ _ ,7^2^ _|_ 



(21) 



Multiplying the integrand by the factor {{q — pY — f^'pi) / {{.1 — vY ~ ™'pi) ^-nd using the 
Feynmann parametrization rule with the change of variable q' = q + p{x — y — 1) — ky, we 
obtain: 

/■^ r f d'^q' q'^ + {kyf + ^pkyiy -x^ + p^ix- yf - ml, 
G{P) = 2i / dx dy — ^ (22) 







(27r)4 {q' + s, 



and 



G{P) 1 r {ky f + 2pky{y - x) +p'^{x -y) 



2 



rrii 



+ — / dx dy'-^ '—^ ^ . (23) 

327r^ In Jn Si + le 



Now, we want to calculate the amplitude for the second diagram in Fig. [T] containing the 
three-vector vertex. The only difference with the previous diagram is the 3V vertex. Thus, 
the amplitudes corresponding to the 3V vertex and the anomalous vertex are, respectively: 

tvv^v, = V,g{{q-p + fc)^6(^')6(^)^e(^^)^ -{p + 2k- qU^^ e^''^^ e^""^^^ 

+ (2(p -q) + k)jPe^^'^>'e^''>} (24) 
tv.j/^P = AG' e-^^\q-p)J^^^p,ef^^ , 

where V3 and A are numerical factors. 

Thus, we can write the amplitude of the diagram (2) in Fig. [l]as 

-it2 = -eGgxC J ^^^^ [q - Pjae'^ 'p^e\ ' e^, 'e^ ^ 
X {(g - p + fc)^e(^')e(^)^e(^)" - (p + 2A; - q^e^^^^'^^^ e^^^^ 
+ (2(p-g) + A:)4^)e(^)^e(^')^} ^ ^ 



^ ' q^ — im?p + ie {q — pY — iriy^ + ie 



{p + k — qY — my + ie 



where G = V3PA. In this process the D*^ is very close to be on-shell with zero three 
momentum. To be consistent with the approach of [26], which is neglecting the three- 
momentum compared to the mass of the vector meson, \q\/mv — 0, e*^^^° ~ 0, we perform 
the sum over polarizations as 

^ |j(Mi^')spatiai ^ ^ij ^25) 

A 

We also can keep the covariant formalism and remember at the end that fi, v' are spatial. The 
way to proceed is very similar to that of the previous diagram. The second term of the 3V^ 
vertex proportional to (p+2/c— g)^ does not contribute, since we have {q—p) aP-yip+^k— q) p = 
qa{p + 2k) pp^e"-^^^ , which applying Lorentz covariance in the integral turns into a term like 
{a'pakjs + h'pi3ka)P'y(-"^'^^ = 0. Therefore, we have two kinds of integrals: 

d^q qa{q-p + k)^, 

{2nY {q^ — m^p + ie){{q — pY — rriy^ + ie){{p + k — q)"^ — my + ie) 
ii^lQau' + hikak^i + CiPak^i + dipyika + eip^'Pa) (26) 
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and 



(27r)^ (g^ — mp + ie){{q — pY — rriy^ + ie){{p + k — qY — my + ie) 

= i{a2gau + hKK + C2PaK + d2Kpu + e2PaPu) ■ (27) 

One can see that only the coefficients proportional to cii, 61, di,a2 and d2 survive. Thus, we 
finally get 

t, = -CeG'gx e"^^^{(aie(f ) + {hk'^ + d,p^)ef^k^)e^J^ + (a^e^ + d2e^;:b'K)ef ^]p4"'^^ , 

(28) 

where now 

G{P) 1 /"^ ^ {kyf + 2pky{y - x) + p'^{x - y f - iTLy^ 



Oi = I dx dy 

bi = [ dx [ dy^^^ ~^ ^ 



16vr2 S2 + ie ^^g) 

di = I dx I dy ^^^ 

167r^ Jo Jo S2 + le 

a2 = —2ai 
d2 = -2di , 

with 

S2 = —my^ + {my^—mp)x + {k'^ — my+m'l,)y+p'^{x — y){l — x+y)+2kyp{x — y) — k'^y'^ . (30) 

In order to evaluate diagrams (3) and (4) in Fig. ([T]), we only have to do the exchanges 
/c <H- p and the amplitudes of diagrams (1) and (2). 

We have 

h = BeG'g'k e^^^^^ae^^^ + dk^{p ■ e^^^))ef^p,ef^^ , (31) 



with 



and 



G{P) 1 f\ f\ {pyf + 2pky{y-x)-m\ 

a = —— + -— / dx / dy — '- (32) 

2 327r^ Jo Jo -ss + te 



d=^,Cdx fdy^^^^, (33) 
Svr^ Jo Jo S3 



where 



S3 = — m,p^ + (mp^ — mp)a; + (p^ + mp — my)y + 2pky{x — y) — p^y"^ , (34) 
for diagram (3), and 

U = -CeG'gx e'^^'''{ia^ei''^ + {c,k^ + e^pnei''^p^)e^/^'^^ + , 

(35) 
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with 



and 



ei 

Cl 
C2 



G{p) 



647r2 



167r .,0 ^0 



dx / (iy 



1 PX 

dx / c/y 
Jo 



S4 + ie 



1 



16vr JO 
-2ai 
-2ci , 



S4 + ze 
y{y-x) 



1 /"X 

(ia; / dy- 
Jo Si + 



2„.2 



S4 = — + (myj — mp)x + {p — rriy + mp)y + 2kyp{x — y) — p y 



for diagram (4). 



(36) 



(37) 



X{P) 



P{q) ^ 4" 
^ I 



J/^{p) 




V(P-q) 



ila) 



\ Pi{q-p) 



p,uj{k) 



X{P) 



P{q) 



V{P - qT 



Vi{q-p) 



(2a) 



J/*(p) 



X{P) 



V{P-q) 




I Piiq-p) 



MP) 



j/<f{p) 



V{P-qT 



p,uj{k) 



Vl{q-p) 



J/^{p) 



(3a) 



(4a) 



FIG. 2: Different types of Feynmann diagrams for the decay of the X(3872) into J/ipp and 

J/ipu. 



This formahsm also allows us to evaluate the amplitudes for the decays X — )■ J/ipp and 
X — )■ J/ipu (Fig. |2]). We can proceed in complete analogy with the radiative decay to 
determine these amplitudes, simply removing the final photon and leaving the vector meson 
in the final state, the p° or the u. Moreover, we must take into account that the and 
the u do not couple to the strange D mesons, so that we have again four different kinds of 
diagrams, but only two channels plus their complex conjugate, that is 16 Feynman diagrams 
to evaluate. Doing this, we can observe that the new amphtudes have the same structure of 
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Diagram P V Pi B 



Diagram P V Vi C 



1 






4 
3^ 


2 











D+D*- 


D+ 


1 

3^2 




D+D*- 


D*+ 


1 

V2 




DtD*- 


Dt 


1 

3v^ 






^ s 


1 


i 






4 
3^2 


2 











D-D*+ 


D- 


]_ 

3^2 




D-D*+ 


D*- 


1 

V2 




D-DI+ 


Ds 


1 

3v^ 




D-D*+ 


^ s 


1 

V2 


3 









4 






4 

3v^ 




D+D*- 


D+ 


1 




D+D*- 




_ 1 
3v^ 




DtD*- 


Dt 


1 




D+Dr 


D*+ 

^ s 


1 

3v^ 


3 









4 






4 
3^2 




D-D*+ 


D- 


1 

V2 




D-D*+ 


D*- 


1 

3^2 




D-D*+ 


Ds 


1 
V2 




D-D*+ 


^ s 


1 

3v^ 



TABLE II: Coefficients B and C of the different diagrams in Fig. [T] 

the previous ones and can be obtained, up to a coefficient, directly with the substitutions 
e g and e^'^-' e^^''^^ For instance, in the case of the diagram (la) of Fig. ^ we have 



= B'gG'gxe'-'^'iael'/^^ + cp^k ■ e^'/^^)ef\,e\''''^ , (38) 
with a and c the same as before 



G{p) ^ 1 J r J i^y)^ + '^p^yi.y - + p^i^ - yY - ^% 



2 327r^ Jo Jo Si + te 

c = /" dx [ dy— — — 
8vr^ Jo Jo si 



(39) 



and B' = PyA of Eqs. ([8]). However, since we are dealing with different vertices, the new 
numerical coefficients, that we call B' and C, are now different and they are written in the 
tables [ml and lYl 
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Diagram P V Fi B' Diagram P V Vi C" 



^0^*0 ^0 1 



D+D*- D+ -\ 



D+D*- D*+ -\ 



^0^*0 ^o_i 



D-D*+ D- \ 



D-D*+ D* 



D+D*- D+ \ 



D+D*- D*+-\ 



D^D*° D° i 



D^D*^ £,*o_i 



D~'D*+ D- -\ 



D-D*+ D*- \ 



TABLE III: Coefficients B' and C of the different diagrams in Fig. |2jin the case of a 

meson in the final state. 



Diagram P V Pi Diagram P V Vi C" 



^0^*0 ^0 1 



D^D*^ D*° \ 



D+D*- D+ \ 



D+D*- D*+ \ 



D^D*^ D*^-\ 



D-D*+ D- 



D-D*+ D*- 



D^D*^ Z)0_i 



D^D*^ D*° \ 



D+D*- D+ -\ 



D+D*- D*+ \ 



D^D*^ ^0 1 



D-D*+ D- \ 



D-D*+ D*- -\ 



TABLE IV: Coefficients B' and C of the different diagrams in Fig. |2j in the case of a 

meson in the final state. 
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III. RESULTS 



Following the procedure in section [IT] we can obtain the total decay amplitude for the 
radiative decay of the X meson and evaluate the correspondent decay width for this channel 
by means of the formula 



\k 



2 



SttMI 



where we sum over the polarizations of the final states and average over the X meson 
polarizations. 

To evaluate the amplitude we use the dimensional regularization for the loop function. 
The subtraction constants used are a = as = —1.27, where the subscript S identifies the 
strange channel, and they are chosen to fit the mass of the X(3872) |26] . 



Applying Eq. (40), we obtain 

T{X Jli)-i) = 149.5 keV . (41) 

In order to make an estimation of the theoretical uncertainty on this quantity we perform 
a suitable variation of the parameters used to compute the total amplitude: the coupling 
G' for the VVP vertex (Eq. (|2])), the axial- vector-pseudoscalar couplings gx for the three 
channels, and the two subtraction constants in the loop function, a and as- 

We allow the constant /, contained in G', to vary, but keeping the relationship Gy = 
f/V2 and replacing My = Mp by Md*. The couplings gx for the neutral and strange 
channels are also varied, independently, by 10%. On the other hand, the variation of the 
coupling for the charged channel is such that the ratio between it and the one for the neutral 
channel is kept constant, in order to preserve the isospin of the X(3872). Then, we let the 
subtraction constants a and as vary between —1.60 and —1.27. This range is motivated by 
the range chosen for /. Indeed, going to higher values of the constant / causes a decrease of 
the potential in the Lippman-Schwinger equation used to evaluate the scattering amplitude 
which determines the position of the resonance. One would need to go to more negative 
values of the subtraction constants a and as in the loop function, which appears in the a 
coefficients, to keep the pole representing the resonance in the same position. The range is 
thus chosen such as to produce an effect in the pole position similar to that induced by the 
change in /. 

We obtain the result 

T{X ^ J/tfj-f) = (117 ± 40) keV . (42) 

We can also evaluate the branching ratios for the decays X — )■ J/ipp and X — )• J/ipu. 
These two decays, if we consider the p and the u with fixed masses, are not allowed because 
of the phase space, but they can occur when their mass distributions are taken into account 
and they are observed in the decays X — > J/ipiTTc and X — )■ J/ipinnT respectively. The two 
and three pions states are produced in the decays of the p and the u. 

Thus, the decay widths, convoluted with the spectral functions, are given by the formula 



1 /•(■'np/„+2rp/„)2 / . 



N 



(mp/ij-2rp/„) 



2 \ TT 



1 



Tx{m)9{mx - mj/^ - rh) 

(43) 
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where 



N 



(mp/^+2rp/„)2 



dm 



Z.2 



TT 



Im 



(44) 



where rx(?Ti) given by Eq. (40), changing nip and m^j by m. 

In Eqs. (43) and (44), rup = 775.49 MeV and = 782.65 MeV are the masses of the 
mesons, Tp = 149.1 MeV and = 8.49 MeK are the on shell widths and 



Qp/u 



(45) 



where q and qp/^ are the off shell and on shell momentum of the mesons in the center of 
mass reference frame: 



4m?. 



-6{m — 2m^ 



Qp/u 



p/w 



, — 4m?. 



(46) 



We call Tx the total decay width of the X into J /ipp or J /ipu to simplify the notation and 
m^r is the pion mass. 
Using Eq. (43) we find 



Vp = 821.9 keV , 
r,, = 1096.6 keV 



(47) 



and when the error anlysis that leads to Eq. (42) is done, the band of values becomes 

Tp = (645 ± 221) keV , 
= (861 ± 294) keV . 



(48) 



With the results of Eq. (47) we can evaluate the ratio 

B{X J/ijmi) Tp 
However, the experiment gives the ratio 

B(X JMvr+Tr-vr") 



(49) 



(50) 



and, to compare our result with this, we must take into account that the uj decays into 
7r+7r^7r° with a branching ratio i?aj,37r = 0.892. 
Hence, our ratio to compare with R^^^ is 



R 



th 



r„ 



X B, 



1.19 



(51) 



well within the experimental error. 
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The result we obtain for the ratio 



^ i"^^' - 0.18 , (52) 



r(X ^ J/^tttt) 



is also compatible with the two values known from the experiment (0.14 ± 0.05) [Ij and 
(0.22 ±0.06) I2j. 



We can also estimate the theoretical errors for the two ratios in Eqs. (51) and (52), 
by evaluating the 7, p and uj decays with the same set of parameters, and varying these 
parameters in the range used to evaluate YiX J /if)'-)): 



= (0.92 ±0.13) 

T(X JM7) , , (53) 

^ (0.17 ±0.02). ^ ^ 



r(X -> J/^Ijtttt) 



The uncertainties in the ratios are smaller than for the absolute values and they are of 
the order of the 15%. 

Finally, we do another exercise removing the D^D*^ — c.c and DgD^^ — c.c and letting 
only the D^D*^ — c.c contribution. The results that we obtain are 

= 0.46 keV 

Tp = 9104.9 keV 

= 368.9 keV , , 

1 54 
i?*^ = 0.04 

r(x ^ j/^7) , 

^ -5.05-10"^ . 



r(x ^ j/v^tttt) 



As we can see, the two ratios that we have to compare with experiment grossy diverge 
from the experimental values and Y p by itself becomes much bigger than the width of the 
X(3872) (Fx < 1.2 MeV). 

In table |V] we compare our results with a variety of results available in the literature using 
different models. It would be interesting to test these models with the new information on 
the experimental ratios to help discriminate among them. 

The ratio of J/V'7 to J/^/'tttt is also evaluated in [9], where the Weinberg compositeness 
condition [T^ is used to determine the couplings but other assumptions are made, and they 
find a range of values from 0.18 to 1.57 depending on the model they consider. We should 
stress that once the X is obtained in our case and, hence, the couplings are determined, there 
is no freedom in our approach, which makes unique predictions, up to the small uncertainties 
tied to the limited range of variation of the parameters. This comparison helps gauge the 
value of the results obtained here. 

We should note that our results are based on the central value of the masses in the PDG 
and a binding of the D^D*^ — c.c oi B = 0.2 MeV . However the error on the mass of the 
X(3872) is as big as that value. This is important to note because the coupling goes as 
according to the Weinberg compositeness condition [13]. When in the future the binding 
can be more accurately determined we can also obtain more accurate values of the absolute 
rates. On the other hand, the values of the ratios will be essentially unaltered. 
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model 



r [keV] 



molecule 
molecule 



cc 



cc 



11(51 
139 [6] 
8|6] 



molecule + cc 
2-+ 



tetraquark 
present work 



cc 



cc 



125 - 250 [8] 
11 - 71 [9] 
2 - 17 [9] 
1.7-2.1 j27j 
45 - 80 do] 

10 - 20 [m 

77 - 158 



TABLE V: Results from previous works for the decay width of the X(3872) into J/'ip'~f, 

using different models. 

IV. CONCLUSIONS 

In this paper we have exploited the picture of the X(3872) as a composite state of 
DD* — cc. dynamically generated by the interaction of the D and D* states. The couplings 
of the state to the different DD* — cc. channels have been calculated before within this 
model and are used here. The coupling for the D^D*^ — cc is similar to the one that would 
be obtained using the compositeness condition of Weinberg, since the state is barely bound 
in the D^D*^ component, but the dynamics of the model produces also couplings for the 
D+D*- ~c.c and DgDg~ — c.c states. Using an extension to SU(4) with an explictit breaking 
of this symmetry of the local hidden gauge approach, used before successfuly in the study 
of related processes, one can determine the widths of the X(3872) to J/ipp, J/ipu and J/ip'j 
and compare with the ratios determined experimentally in recent works. We find a very good 
agreement with the experimental results. The absolute numbers obtained for the different 
widths are also reasonable and their sum within errors, (1.6 ±0.6) MeV, is compatible with 
the recent total X(3972) upper limit of the width, F = 1.2 MeV. 

We have also conducted a test neglecting the charged and strange components of the 
wave function and thus having only the D^D^* — cc component. We obtain ratios in great 
disagreement with experiment and an absolute value for the X(3872) partial width into J/ ipp 
which largely exceeds the experimental upper bound for the total width of the X(3872). This 
exercise confirms the relevance of the charged channels and the approximate 1=0 character 
of this resonance, contrary to the appealing picture of the D^D^* — cc state, which has a 
larger probability to be found because it is less bound than the charged components. We 
showed that in these processes it is the wave function at the origin what matters, or more 
concretely the couplings, which are related to it, and not the probability. 
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